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The hnear stabihty of electrically driven flow of liquid metal in circular channel in the presence of 
vertical magnetic field is studied. It is shown that the instability threshold of such flow is determined 
by magnetorotational instability of non-axisymmetric modes (m 7^ 0) and does not depend on the 
type of the fluid if magnetic Prandtl number is small Pr <C 1. Our numerical results are found 
to be in a good agreement with available experimental data from Grenoble High Magnetic Field 
Laboratory, France [P. Moresco and T. Alboussiere, J. Fluid Mech. 504, 167 (2004)]. 



The theoretical and experimental study of magnetoro- 
tational instability (MRI) has attracted much attention 
in recent years. MRI was originally predicted by Ve- 
likhov in 1959, but its intensive study began only in 
1991 when it was rediscovered in astrophysical context 
by Balbus and Hawley At present time subject of 
MRI is one of the most important developments in mag- 
netohydrodynamic (MHD) theory with far reaching con- 
sequences for a varietyof astrophysical phenomena such 
as accretion disks magnetic reconnection Q and 

dynamo Q. 

One of the topics of great current interest is experimen- 
tal verification of MRI. Several experiments have been 
initiated to investigate MRI in laboratory 
[HI . [T^ by studying the stability of conducting fluid (liq- 
uid metal) rotating in transverse magnetic field. Despite 
these attempts, the MRI has never been clearly detected 
in laboratory and any progress in this direction is ex- 
tremely important. Current status of experimental stud- 

ll[l4|. 



ies of MRI is described in 



Two different mechanisms of the fluid rotation have 
been proposed for MRI experiments so far: mechanical 
drive by virtue of viscous drag force acting on the fluid 
between moving surfaces (Couette flow) and electrical 
drive by Ampere force arising when the electric current 
is passed through the fluid in transverse magnetic field 
(electrically driven flow). In most existing MRI experi- 
ments a Couette flow is used either in cylindrical 0, d, Q 
or spherical geometry [l3|. The main difficulty of the 
cylindrical Couette flow is the presence of the stationary 
end-caps that affect the entire equilibrium flow making it 
different from the idealized infinite-cylinder angular ve- 
locity profile ri(r) = a + b/r"^ so the conditions for MRI 
may not be met. Their influence can be reduced either by 
employing the differentially rotating end-caps 14, or 
by using sufficiently long cylinders. In the latter case, 
experimental observations of MRI have been reported 
[^, [1^ . MRI has also been observed in the experiment 
with rotating spheres [l3| though the background flow 
was already fully turbulent without any magnetic field, 
indicating that MRI is not the only possible instability 
in this geometry. 




FIG. 1: Sketch of electrically driven flow. Some details are 
shown as in Fig. 1 in Ref. [l9l ]. 



Another way to rotate conducting fluid in circular 
channel with axial magnetic field is to apply radial elec- 
tric current (Fig. [T|). In this case the equilibrium flow 
forms so-called Hartmann layers near the top and bot- 
tom walls and parallel boundary layers near the side 
walls. The widths of these layers scale with Hartmann 
number Ha (see Eq. ^ for definition) as 0(Ha~^) and 
0(Ha~^^^) respectively, so they become negligible at high 
values of magnetic field [l3| • Such flow has the angular 
velocity profile i^(r) oc almost entirely in the cross- 
section of the channel [l8|. This profile is stable in hy- 
drodynamics according to Raleigh's criterion, but it can 
be destabilized in the presence of magnetic field since it 
satisfies the necessary condition for MRI dil^{r)/dr < 
[l|. MRI experiment based on electrically driven flow of 
liquid sodium has been proposed in the Russian Research 
Center "Kurchatov Institute" [ll| and built in Obninsk 
(Russia). A similar configuration using plasma instead 
of liquid metal has been developed in Los Alamos [l^ . 
At present time experimental data are not available from 
those MRI experiments. 

The electrically driven flow of liquid metal (mercury) 
in circular channel was also used in Grenoble High Mag- 
netic Field Laboratory, France to investigate the stability 
of the Hartmann layers . In this experiment a well- 
marked transition to turbulence was found when the ratio 
of Reynolds number to Hartmann number R = Re/Ha 
exceeded critical value R* « 380 (star denotes the quan- 
tities taken from [l^l ) , this value of R* was valid also for 
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TABLE I: Parameters of Grenoble and Obninsk experiments. 



Ri 


(cm) 


i?2 (cm) 


L (cm) 


Fluid 


Pr 


Grenoble 


4 


5 


0.5 


Hg 


1.49 • 10- ^ 


Obninsk 


3 


15 


3 


Na 


8.8-10-6 



inverse process of laminarization without a visible hys- 
teresis and for a wide range of intensities of the magnetic 
field (Ha = lOO-MTOO). The critical value R* obtained in 
Grenoble experiment is two orders of magnitude smaller 
than K^'^ w 50,000 predicted by the linear stability the- 
ory for Hartmann layers [20| . A possible explanation to 
this paradox is discussed in papers [21, 22|, where a com- 
plicated two-step transition scenario to turbulence in the 
Hartmann layers is assumed. In these papers the radial 
dependence of the velocity profile in the main part of the 
flow is completely disregarded. In the present Letter, it 
is suggested however that in this configuration the mag- 
netized flow can be destabilized via the MRI mechanism. 

We study the linear stability of electrically driven flow 
of liquid metal in the geometries relevant to the Grenoble 
and Obninsk MRI experiments (Table [I| and show that 
the instability thresholds in these experiments are deter- 
mined by MRI of global non-axisymmetric modes with 
large azimuthal mode numbers m. In our analysis we 
neglect all boundary layers and assume that the fluid ro- 
tates in the uniform axial magnetic field Bq = BqGz with 
equilibrium velocity Vq = rf2(r)e(p where ^{r) = Mq/v^ 
throughout the cross-section of the channel. The equi- 
librium angular momentum of the velocity A/q is deter- 
mined by the total electric current /q passing through the 
channel llSl: 



where p, a and v are fluid density, electric conductivity 
and kinematic viscosity respectively. 

A number of similar studies was performed by different 
researches 0, S, [H, [l^l , who considered the stability of 
more general idealized Couette flow ri(r) = a + b / r"^ . The 
major part of these studies is restricted to the case of ax- 
isymmetric modes with to = 0, S] , which are known to 
have the largest growth rate of MRI. Numerical analysis 
of non-axisymmetric modes with small to (23l . [2^ shows 
that they can have lower instability threshold (lower crit- 
ical values of Reynolds number Re), so they might be eas- 
ier to excite in real experiment. For electrically driven 
flow it was found recently [2^ that in ideal MHD the 
MRI threshold decreases with to as Re oc 1/m; therefore 
the most dangerous modes in such flow are those with 
larger to (though they have smaller instability growth 
rate). Our present results suggest that overall stability 
of the flow in a finite height circular channel is deter- 



mined by dissipative perturbations with large azimuthal 
numbers to. 

We use dissipative incompressible MHD equations lin- 
earized about the equilibrium state and represent all per- 
turbations in the form /(r, z) exp(7< -|- irmp) in cylindri- 
cal system of coordinates {r, tp, z}. It should be stressed 
here that we consider the channel of the finite height - 
this is a substantial difference of our stability analysis 
from previous studies. For computational convenience 
we introduce dimensionless quantities, taking as a unit 
of length a half-height of the channel L and as a unit of 
time the characteristic viscous time I? jv. The perturba- 
tions of velocity i5V and magnetic field (SB can be written 
in terms of dimensionless vectors v and h as 
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where 



Ha = 



' — , Re = , 

pv v 
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arc Hartmann, Reynolds and magnetic Prandtl numbers 
respectively. Introducing vortex u; = V x v we arrive at 
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w^e^ + HaV X h^,(2) 
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hr&^p + Hav'^, (3) 



where prime denotes the derivative with respect to z. 
Note that r and components of these equations are 
enough to find the full solution (2-components can be 
deduced from the conditions V • u; = and V • h = 0). 

The proper boundary conditions should be specified 
prior to solving the system ([2]), In the flow of viscous 
fluid all velocity components vanish at the rigid walls, i. e. 



-0, 



=1 = 0, 



(4) 



The boundary conditions for magnetic field depend on 
the conductivity of the walls. In both Grenoble and Ob- 
ninsk experiments the side walls of the channel can be 
considered as perfect conducting. At the surface of the 
perfect conductor the time- varying normal component of 
magnetic field as well as the tangential components of 
electric current should be zero. This means 



hr\r—ri,r2 — 



d{rh^) 



dr 



0. 



(5) 



The other two walls (Hartmann walls) are electrical in- 
sulators. For simplicity we assume that the perturbed 
components of the field are zero at these walls, i. e. 

/ir|z=±l = 0, h^\z = ±i=0, (6) 

which is consistent with the absence of the normal com- 
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FIG. 2: Calculated marginal stability curves for geometry 
of Grenoble experiment. Azimuthal mode numbers m are 
shown. Difference between Hg and Na is within the line width. 
Straight line corresponds to the transition to turbulence ob- 
served in experiment Re = 1, 700Ha. 



poncnt of the current at the surface of insulator. 

Equations ([2]), ^ with boundary conditions (01)-® 
constitute an eigenvalue problem with 7 being an un- 
known eigenvalue. A solution to this problem was sought 
by expanding functions Vr, v^p, hr and in terms of ei- 
ther odd or even polynomials in z (up to Nz = 8), and 
by discretization of the system ([21), © in terms of finite 
differences in r-direction (up to Nr = 50). Then this 
system was reduced to a large {ANrN^ x ANrN^) ma- 
trix eigenvalue problem which was solved using standard 
numerical methods of MATLAB. 

The following numerical procedure was performed for 
two values of magnetic Prandtl number corresponding 
to mercury (Hg) and liquid sodium (Na) and for the 
geometries of both Grenoble and Obninsk experiments 
(Tabic Taking a value of azimuthal number from the 
array m = -f- 200, we scanned through a range of val- 
ues of Re and Ha, finding the maximal growth rate for 
given parameters. For each value of Ha we determined 
the value of Re that yields a marginal stability, i. c. cor- 
responds to the zero maximal growth rate Re(7) = 0. 
Thus we obtained the marginal stability curves at the 
plane Re - Ha (Fig.[a[3l). 

One can see from Fig. [2] and [3] that the instability 
threshold is determined by different azimuthal modes at 
different Hartmann numbers Ha. For larger Ha the corre- 
sponding m is larger. In fact, a simple scaling law can be 
obtained for large m: Ha oc m^. The stability curves of 
axisymmetric modes with m = are not shown because 
they are situated at much higher values of Re. There- 
fore, the non-axisymmetric modes play the decisive role 
for excitation of MRI in electrically driven flow. 

It is worth noting here that the marginal stability 
curves practically do not depend on the magnetic Prandtl 
number Pr of the fluid; they are determined only by 
the geometry of the channel. This is true for non- 
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FIG. 3: Calculated marginal stability curves for geometry of 
Obninsk experiment. Azimuthal mode numbers m are shown. 
Difference between Hg and Na is within the line width. 



axisymmetric modes in the limit Pr <^ 1 which is com- 
mon for most liquid metals. In this case a good approx- 
imation for instability threshold can be achieved by ne- 
glecting all the terms containing Pr in the Eq. ([3|) and 
reducing the system ([2]), ([3]) to one vector equation with 
hydrodynamical variables. Such approach is not applica- 
ble for axisymmetric modes in axial magnetic field [26| . 

In ideal MHD, the MRI threshold to a significant de- 
gree is affected by singularities inherent in eigen-value 
problem 



27|. In incompressible limit, the MRI of 



non-axisymmetric modes is associated with one set of 
these singularities - the so-called Alfven resonances jisf . 
Within the frame of MHD model considered here, the 
Alfven resonances arc removed by dissipative effects (re- 
sistivity and viscosity). An important consequence of 
this is that the dissipative stability threshold for non- 
axisymmetric modes appears to be lower than the ideal 
one, i. e. the dissipation (mainly the resistivity) has 
destabilizing effect on the ideal modes. 

Our main observation is that the MRI threshold of the 
electrically driven flow is formed by the envelope of all 
marginal stability curves corresponding to modes with 
different azimuthal numbers m. The shape of the enve- 
lope depends on the particular geometry of the channel 
(see Fig. [2] and [3|) . In the Fig. [21 the envelope is close 
to the form Re oc Ha when Ha ^ 300. It means that 
instability should be excited in this geometry if the ra- 
tio R ~ Re/Ha exceeds some critical value Rc. As men- 
tioned above this effect was actually detected experimen- 
tally 

For comparison of our numerical results with the ex- 
perimental data we need to calculate Rc taking into ac- 
count our definition of Reynolds number ([T]). The rela- 
tion between Reynolds number Re* from Ref. [l^ and 
Re used in our calculations is 



Re* 



dv„ 

V 
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turbulence observed in the experiment [l9[ is associated 
with non-axisymmetric MRI in electrically driven flow. 

This work is supported in part by NSERC Canada. 



FIG. 4: Calculated contours of marginally stable eigenfunc- 
tion in cross-section of Grenoble channel at m = 100, 
Ha = 500, Re = 9 ■ 10^. Magnitude is shown in arbitrary 
units. 



where d ~ 2L is the channel height and Vm is the mean 
velocity in the equilibrium flow. Our definition of Ha is 
the same as definition of Ha from Ref. [l^ used in the 
figures (though it is not the same as stated in Eq. (2.4) 
in Ref. [l9[). Thus we obtain 



" Ha ~ 0.22 



1700. 



The line associated with this experimental value of Rc 
is also plotted in Fig. [2l As one can see, the calculated 
threshold of MRI is in a good agreement with the exper- 
imental results. 

A natural question arises: what do we observe in re- 
ality - MRI or Hartmann layer instability? A strong 
argument for MRI is the comparison of respective linear 
instability thresholds: MRI threshold found in our cal- 
culations corresponds to that measured in the Grenoble 
experiment and two orders of magnitude smaller than 
instability threshold of Hartmann layers [2^. Also it is 
evident that in the presence of global (affecting the entire 
flow) robust linear instability such as MRI the nonlinear 
effects in Hartmann layers |2l|, are unlikely to play the 
major role in destabilizing the flow. The global character 
of MRI is illustrated in Fig. 2] where a typical marginally 
stable eigenfunction v^p is shown. 

In conclusion, we have shown that the instability 
threshold in electrically driven flow is determined by MRI 
of global non-axisymmetric modes with large azimuthal 
mode numbers m. This threshold does not depend on the 
type of conducting fluid as long as the magnetic Prandtl 
number of the fluid is small Pr ^ 1. MRI threshold 
calculated for the geometry of the Grenoble experiment 
agrees well with the critical ratio Rc found in this ex- 
periment. These results suggest that the transition to 
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